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ABSTRACT 


Recent .interest in asymmetrically swept, or oblique, wings has raised 
fundamental questions about the aoroelastic stability characteristics of 
such wings. This paper presents two formulations of the oblique wing 
flutter problem; one formulation allows only simple wing bending deform- 
ations and rigid body roll as degrees of freedom, while the second 
formulation includes a more complex bonding-torsional deformation together 
with the roll freedom. Flutter is found to occur in two basic modes. The 
first mode is associated with wing bending-aircraft roll coupling and occurs 

at low values of reduced frequency. The second instability mode closely 
♦ 

resembles a classical bending* 'torsion wing flutter event. This latter mode 
occurs at much higher reduced frequencies than the first. The occurrence 
of the bending-roll coupling mode is shown to lead , to lower flutter speeds 
while the bending-torsion mode is associated with higher flutter speeds. 

The ratio of the wing mass moment of inertia in roll to the fuselage roll 
moment of inertia is found to be a major factor in the determination of which 
of the two instabilities is critical. 
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Nomenclature 


[a] “.flexibility matrix for clamped fuselage wing 

b =* wing reference semi-chord 

r 

c « wing chord measured perpendicular to elastic axis (Fig* 1) 
“ wing chord measured parallel to the free stream direction 

r 

o- « 2-dirtcnsional sectional lift-curve slope 

g = structural damping parameter 

i 

“ Wing roll moment of inertia at zero sweep 
k = reduced frequency, mc/2V^ or toc/2V 
R “ (Iq/I^) cos^A . 
s = wing semi-span 

V = airspeed 

flutter speed . ; 

“ airspeed normal to swept axis, “ VcosA 
A = sweep angle 

p -air density 

0 ) = frequency of oscillation 


t 
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Introduction 

The recent interest in the use of an asymmetricany swept, high- 
aspect-ratio wing to achieve high 1 i ft- to-drag ratios has generated in- 
terest in the aeroelastic stability characteristics of such a configuration. 
However, the undesirable static aeroelastic divergence characteristics of 
symmetrically swept forward wings has prompted some caution on the part 
of structural engineers towards the asymmetrical wing. As a result, con- 
siderable discussion of the merits of such a design and the potential 
weight penalties which might be incurred has occurred. Oones and Nisbet 
•{Ref, 1) have presented data which tend to allay some misgivings about 
the aeroelastic stability of asymmetrically swept or oblique wings. Pro- 
minent among their findings is the discovery that the inclusion of the 
rigid-body roll degree of freedom of the aircraft appears to have a stabi- 
lizing effect oh the aeroelastic stability of the wing, when compared to 
the stability of a similar, but clamped, wing. Their analytical results 
were obtained through the use of quasi -static aerodynamic th-^ory to re- 
present the perturbation lift forces generated by the harmonic motion of 
their idealized flexible model. 

This study seeks to explore, in somewhat more detail than Ref. (T),' 
the flutter behavior of asymmetrically swept or oblique wings; to accom- 
plish this task the results of two studies are presented. The first 
study examines the flutter behavior of an idealized finite span, uniform- 
property wing in .incompressible flow swept asymmetrically at various angles 
to the flow. For this portion of the study, quasi -steady aerodynamic strip 
theory will be employed in the equations of motion; the Galerkin method 
will be used to solve these equations. 
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The second portion of tudy entails the use of a more sophisti- 
cated approach to the solution of the oblique wing flutter problem. This 
approach uses a finite-element, unsteady aerodynamic representation 
together with a multi -degree-of-freedom structural model to examine more 
closely and more accurately the flutter behavior of variable ixfiform 
wings, In all cases, the flow is assumed to be incompressible, 

From these studies, it will be shown that, at moderate sweep angles, 
the flutter speed of the wing may be lowered when compared with the flut- 
ter speed of the wing at zero sweep. In addition, the shape of the wing 
planform and the spanwise distribution of stiffness and weight will have 
a significant effect on the relation between flutter speed and sweep 
angle. 

Discussion 

The first part of this study is concerned with the aeroelastic anal- 
ysis of a simplified oblique wing model, shown in Fig. 1. The impetus 
for such a study stems from the desirability of assessing the behavior of 
the flutter speed of the wing as it is asymmetrically swept. This model 
represents a wing of uniform structural and aerodynamic properties, asym- 
metrically swept at an angle a to the flow. This high-aspect-ratio wing 
is idealized as a beam with a straight elastic axis, free to roll about 
an axis parallel to the flow. It is assumed that mass is distributed 
along this roll axis such that a mass moment of inertia, I^, simulating 
the roll moment of inertia of the fuselage, appears concentrated there. 

To examine the aeroelastic stability of this model, assume that it 
is caused to undergo small oscillations about a "wings-level" static 
equilibrium position. The stability of the subsequent motion can be 
determined by an examination of the character of this free vibration. 
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The structural behavior of this wing can be modelled through the use of 
conventional Euler-Bernoulli beam theory. It is further assumed that the 
wing has no torsional flexibility so that only bending flexibility is 
important. The limits to the validity of this latter assumption will be 
discussed later in this paper. 

In Ref. 2 , Barmby> et al . discuss the flutter analysis of 
symmetrically swept wings through the use of aerodynamic strip theory and 
the Theouorsen ' functions. The present study neglects all the noncircu- 
latory aerodynamic terms in Ref. 2, but retains two of the circulatory 
terms. In addition, the free vibratory motion is assumed to take place 
at a value of reduced frequency k = iiic/2V,^ which is so small that the flow 
is quasi -steady. The circulatory aerodynamic terms retained are: a term 
which corresponds to the familiar damping-in-rolljand a term which 
arises from the angle* of attack generated by bending deformations of a 
swept wing. 

The assumptions about the behavior of this idealized model under- 
going small oscillations in’ the airstream lead to the following 
differential equation of motion for the elastic wing. 

m ^ + El ^ + (qcc^^cos^A), ||- tanA (1 ) 





2 

^ qcc. cos A 
3t “ 



- mpycosA = 0 

where m = wing mass per unit length along the y-axis 

El = bending stiffness of cross-section perpendicular to y-axis 
q = frees.tream dynamic pressure 

W = wing deformation due to elastic deformation, positive upward 


t = time 

p - roll rate in radians per unit time 
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NoncHmensionalization of Eq. (1) yields the following equation. 


w + A ^ , aw ^ f n 

IT ) ^ ^ ^ On ItanA jlv 


( 2 ) 


^wl\ 


4-) (pcosA) n = 0 


^sin /W I \El 
whore (’) = differentiation with respect to time 
w == W/L 
n- y/l 

X = qccj^^L sinAcosA/EI 

The requirement that the sum of all roll, moments generated by wing oscil- 
latory motion be equal to zero results in the additional equation: 


[I.J. + IqCos^a ) p = ( qccL^L^cos^A^ ^ 


ndn ) tanA 


+ I ml'^cosA 


1 .. ■ 

wndn 


p o / \ 

. 2 qcc. L cos A / pL \ 

”3 V / 

t 

qcc. L^cos^A , 

y.iwndn 

If we let p cosA then Eq. (3) may be written as 


+ 3 Iv 




(ysinA) 


/-I s -!(«/; 


1 .. 
Wndn 




Wndn 


where Ij = t IqCos^A 


2 2 3 2 

I, = Lcos A = l-mL^cos^A 
W 0 3 

2 3 

Y = qccj^ L cos A/Ij 


(3) 


(4) 
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To solve Eqs. (2) and (4), the time dependency is eliminated by recogni- 
tion that. the functions w{njt) and <ii(t} are separable sfich that 

w(n,t) = f(n)e’"^ (5a) 

Next, Eq. (2) is separated into two parts, one valid in the region -1 n 
5 0, the other valid in the region 0 < n <1, Finally the resulting set 
of equations is solved approximately through use of GalerKin's method. A 
simple polynomial to use fur such a solution is that shape obtained for 
uniform loading of a cantilever beam. In this case the function f(n) is 
approximated as: 

where a and b are unknown constants. The Galerkin method leads to a set 
of three homogeneous algebraic equations, represented in matrix form as: 


r -1 ( a 

0 


dii b 

= 0 

(7) 

L 1 $ 

0 



The coefficients d. . are given in, the Appendix to this paper. 

1 J 

It is found that, in the absence of the roll freedom, the first 
natural frequency of vibration of the clamped wing, in vacuo , is pre- 
dicted by the Galerkin method to be 

“o = 3.530 Jir (8) 

This compares with the exact solution (Ref.. 3) 

(0^ = 3.518 /eT" ■ (9) 

1 mL^ 

For the clamped wing, it is found that the sweptforward wing undergoes 
static divergence whenr is zero. This occurs at a value of X equal to 


(6n^ - 4n^ + n^} 0 - n ~ 1 

j (6n^ + 4n^ + n^) - n “ 0 
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6.40. The exact f?olution cjlves a value of A for static divergence of 
6.33 (Ref. 4). 

If all the system parameters, such as El, A and V, are substituted 
into the expressions for d^.j, then the determinant of the matrixjd^.J, 
written as A(d..), can be used to find r through the relation: 

^ J 

A(d.j.) = 0 (10) 

With reference , to Eqs. (5a, b), it is seen that if r is found to be a real 
number, then motion is aperiodic. A positive real value of r indicates 
aperiodic instability or static divergence. On the other hand, if r is 
found to be a complex number, motion is harmonic. If r = « + ioj then the 
motion is periodic with frequency u. For negative values of «, the motion 
decays, but for positive values of « it grows with time. This latter 
situation corresponds to the dynamic aeroelastic instability commonly 
referred to as wing flutter. At the value r = iw, the system undergoes 
undamped oscillation and is said to be in neutral equilibrium. For a 
given set of system parameters, the airspeed V at which this occurs is 
termed the flutter speed, Vp. 

The way in which the problem is presently formulated allows the selection 
of one system parameter as the unknown in Eq. (10). The magnitude of 
the complex number r is of no interest, but rather the value of velocity 
at which neutral stability occurs. For this reason, it is found to be 
advantageous to let r = iw in the expressions for d.. and to express these 

1 J 

coefficients in terms of w^ and the parameter B, defined as 

3 = A/Aj.y (11a) 

where A^^.^ = 32/5 = 6.40 _ (11b) 
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Given the system physical parameters, the determinant in Eq. (10) 
may be expressed in term'; of the independent variable (1. Collecting 
terms, the determinant is found to have a real part and an imaginary part 
given respectively by the expressions 




t[i- 


0 + BR + 41 ,|;d 
25 6ir . 


= 0 


and 




where 


and 


+ | (1 + 6^/25) = 0 



(13a) 

(13b) 


The selection of a value of 3 which yields identical values of the 
ratio w/Wq in Eqs, (12a, b) completes the solution. With this value of 3» 
the flutter velocity then may be obtained. 

The above solution procedure v/as implemented for a small model v/ing 
constructed of aluminum sheet with a constant thickness of 0.064 inches. 
The wing properties were taken to be: 

Material density = 0.101 Ib^^^/in? 
c = 4 in. • L = 20 in. Cj^^ “ 2it 

y = 3 El = 874.0 Ib^-in. 

Using a sea level air density value, Eqs. (12a, b) were solved numeri- 
cally using a Wewton's Method trial and error solution technique. The 
results of this analysis are shown in Fig. 2. 
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From Fig. 2 , it is seen that the flutter speed decreases as the wing 
is swept, For small values of A, the value of Vp greatly exceeds that of 
the clamped-wing static divergence speed, Vq, However, as A increases, 
the critical speeds Vp and draw closer together; at A = 90® they will 
coincide. As suggested by Jones and Nisbet in Ref. 1, the moment of 
inertia ratio lyi-j- plays a significant role in the flutter analysis ofthis 
asymmetric wing. From the expression for given below Eq. (4), it is 
seen that this ratio tends to zero as A approaches 90®. It has been sug- 
gested that this mass moment of inertia ratio should be as large as 
possible to improve flutter performance. The results in Fig. 2 support 
this C'b^ervation. 

Since one of the original assumptions of the present analysis was 
that the flutter instability occurs at relatively small values of reduced 
frequency k,it is worthwhile to note the values of reduced frequency for 
which' the instabilities in Fig, 2 occur, These numbers are listed in 
Table I. Although these reduced frequencies are reasonably small, the 
accuracy of these results is probably degraded somewhat by the quasi- 
steady flow assumption. 

The model just analysed is similar to, but not identical to, a 
series of models used by Papadales (Ref. 5) in wind tunnel experiments at 
Virginia Polytechnic Institute. Those experiments had as their primary 
objective the study of the static aeroelastic characteristics of clamped 
oblique wings. However, when those tests were completed, simple flutter 
tests were conducted on roll -free models. Although no attempt to take 
accurate data was made during these flutter demonstration tests, the 
velocity magnitudes shown in Fig. 2 correspond to the order of magnitude 
of the velocities observed in these demonstrations. In addition, for 
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sweep angles greater than 15-20°, the primary mode of instability was 
observed to be primarily a fundamental symmetrical bending mode coupled 
with the asymmetrical rigid body roll oscillation. 

For sweep angles less than about 20°, the tests described in Ref. 5 
found a flutter mode which resembled a more conventional "fixed-root" 
bending-torsion instability. These observations, together with the desire 
to obtain a more accurate versatile analysis model, suggested the application 
of a more sophisticated analysis technique to the oblique wing flutter 
problem. , It is to this analysis that attention is now turned. 

Conventional flutter analysis of realistic aircraft employs assumed 
structural deflections or mode snapes together with generalized coordi- 
nates assigned to' these mode shapes. An excellent discussion of modal 
and npn-modal matrix methods of flutter analysis iS’ given by Rodden in 
Ref. 6. In addition. Ref. 6 presents a succinct discussion of how to 
include free-free boundary conditions into the conventional restrained or 
clamped model. This latter discussion follows the development given in 
Ref. 7, but is more general. The highlights of Ref. 6 are reviewed here. 

To analyze tne flutter behavior of a planform such as that shown in . 

Fig. 4, it is necessary that the following items be taken into account: 
the distributed mass of the wing along the span; the variable bending and 
torsional stiffness along the span; and the unsteady, three-dimensional 
aerodynamic forces and moments associated with deformations caused by 
wing oscillations. With the assumption of simple harmonic motion at fre- 
quency to, the classical matrix equation for flutter analysis, before the 
inclusion of assumed modes, is given by (Ref. 6): 

” 1 =(t^) H [W ^ M]H 
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In this equation the static flexibility matrix [a] has been divided by 
the factor (1 + ig) to account for the structural damping necessary to 
sustain simple harmonic motion. The elements of the vector {h} are 
actual elastic deflections and rotations at control points on the wing. 

The mass matrix [M] and aerodynamic influence coefficient matrix (Cj^] are 
both multiplied by the frequency squared. The elements of [Cj^l are 
complex numbers and functions of Mach number and the local control point 
reduced frequency, k = wb/V, where b is the local semi-chord. With the 
formulation in Eq. 14, the unsteady aerodynamic forces enter into the problem 

mathematically, as complex masses. 

The idealization of the wing structure as an assemblage of beams, 
each with a straight elastic axis, permits the use of conventional finite 
element structural analysis methods to describe the wing stiffness and 
flexibility. The reader is referred to Refs. 8 and 9 for discussions of 
this metho'' Similarly, the mass matrix may be formulated from finite 
element methods. The mass matrix must account for the fact that the wing 
shear centers may be offset from the wing chordwise location of the cen- 
ters of mass. Finally, to model the three-dimensional aerodynamic forces 
and moments, a doublet-lattice method (Ref. 10) generate 

aerodynamic influence coefficients. To expedite these calculations, an 

existing computer program (Ref. 11) was used. 

A computer'-program was written to calculate the matrices'in Eq. 14. 

The' free vibration modes for the clamped system are then used to condense ■ 
the matrix equations. The free-free boundary conditions are then intro- 
duced to "free" the clamped system described in Eq, 14; this 
allows- rigid body rolT freedom. Once these matrices have been formed, 
the eigenvalues and eigenvectors may be found. Since the aerodynamic 
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influence coefficients are functions of reduced frequency k and Mach number 
(in these, studies, Mach number is zero), a set of eigenvalues and eigenvectors 
corresponding to each value of k is generated. The familiar V-g method (Ref. 
3, pp. 565-568) is then used to find the value of velocity at which neutral 
stability occurs. 

To assess the effect of torsion and unsteady aerodynamics on the 
flutter analysis of the oblique wing, the uniform property aluminum wing was 
again analyzed. The wing is considered to have the same structural pro- 
perties as before,' but,. in the present example, GJ is taken to be equal 
to 1346 lb-in. It should be noted that the flat, sheet-aluminum wing has 
a ratio of first bending to first torsion which is slightly higher than 
that common to conventional aircraft. 

The analysis of the constant property wing, including roll freedom 
and torsional flexibility and employing the doublet-lattice method was 
conducted with a sixty degree-of-freedom model. These sixty degrees of 
freedom were obtained by considering ten control points on each wing; 
each control point has pitch, plunge and bending rotation elastic degrees 
of freedom. This model was subsequently reduced to a twenty degree-of- 
freedom model by using the first twenty natural modes of the system. 

The results of this flutter analysis are displayed in Fig. 3 as 
ratios of the instability velocity (either "Tutter or divergence) to. the 
velocity at whicn wing torsional divergence occurs at zero sweep; this 
latter velocity 'is denoted as Vpg. 

. In Fig. 3 the behavior of the wing when the fuselage is clamped is 
shown as the curve labelled Vg/VgQ. With the fuselage clamped, instabi- 
lity is found to occur at a reduced frequency k = 0; this is divergence. 

When roF freedom is allowed, and when I^/I^ = 3, a dynamic instability 
appears;, this is flutter and is shown as the curve Vp/V^g-. Unlike the 
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previous results, the flutter speed does not tend to infinity as A tends 
to zero. This analysis reveals that the wing has two possible modes of 
flutter instability. The first. type of instability is characterized by a 
classical bending-torsion oscillation of the wing with its root fixed or 
• damped. The second mode of instability is one which involves bending- 
torsion deformation coupled with a significant amount of rigid body roll. 
Depending upon the sweep angle, one of these modes will occur before the 
other as airspeed. is increased. 

For the example shown in Fig. 3, the bending-torsion-roll mode does 
not become Critical until the sweep angle is near 15°. The cusp in the 
Vp vs. A curve indicates that the mode of instability changes at this point. 
Also, the flutter frequency will change discontinuously at this point. 

■ Gaukroger (Ref. 12) has noted a similar phenomenon associated with 
the pitching degree of freedom on bilaterally symmetric aircraft. A 
phenomenon termed "body- freedom flutter" is found to occur if the aircraft 
pitching moment of inertia is small enough. For large values of pitch 
inertia, a type of wing flutter resembling a "wing-clamped" instability 
occurs. This latter type of flutter is termed "fixed-root flutter." The 
flutter frequency at which body-freedom flutter occurs is significantly lower 
than that at which fixed- root flutter occurs. It is also shown that for 
bilaterally symmetric aircraft anti symmetrical body-freedom flutter involving 
rigid-body roll, is theoretically possible,, but only for cases where the 
fuselage roll moment of inertia is extremely small or. negative in comparison 
to that of the wing. These results are contrasted with the present results 
in which body-freedom flutter occurs if the fuselage moment of inertia is 
large when compared to that of the wing. 
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As a furtb" ■■ illustration of the flutter behavior of oblique wings, ‘ 
a non-uniform wing planform, constructed of the same material as the uni- 
form property wing, was analyzed. This wing {Fig. 4) has a modified 
elliptical planform. In this case, the wing-fuselage combination has a 
•roll moment of inertia ratio ~ 11.69. Fig. 5 shows the stability 

behavior of the clamped and roll-free wings. While the decrease in divergence 
speed with increasing A for the clamped wing shown in Fig. 4 resembles 
that of the uniform property wing, the behavior of the flutter speed for 
the nonuniform wing is much different. Once again, for large sweep angles, 
the decrease of flutter speed with sweep angle is seen in Fig. 5; however, 
the roll -free flutter speeds and clamped divergence speeds are more widely • 
separated in Fig. 5 than in Fig. 3. . 

■ ■ To assess the obvious importance of the fuselage roll moment of inertia, 
the aeroelastic stability of the nonuniform property wing shown in Fig, 4 is 
again studied. However, the roll moment of inertia of the fuselage is now 
increased by a factor of two. The results of this study are presented in 
Fig. 6 and are compared to those previously obtained using the smaller 
fuselage roll moment of inertia. Once again,, the results are displayed as 
ratios of flutter speed to clamped divergence speed at zero sweep angle.- 
The effect of increasing the fuselage roll moment of inertia is 
clearly seen in Fig. 6. The flutter speeds for both moment of inertia 
ratios are seen- to be nearly identical until about 15 degrees- of sweep. 

Near this point, the flutter mode for the Iq/I^ = 5.85 wing changes from 
the fixed-root type to the body-freedom type. This is seen to depress the 
flutter speed as A increases. 
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As a final example, consider the uniform property dluminum wing. This 
wing has been previously analyzed using a quasi-steady strip theory model 
with elastic bending degrees of freedom and roll coupling. It has also 
been analyzed with a bending-torsion model which used the doublet-lattice 
aerodynamic loads. For the present example, the value of the torsional 
stiffness GO is chosen to be 10 times that of the example whose results 
were presented in Fig. 3. .The results of the present study are shown in 
Fig. 7, where they are compared with those presented in Fig. 3. In Fig. 7, 
the designation "Wing 2" refers to a uniform property wing with properties 
identical to those of "Wing 1" except that the wing sectional torsional 
stiffness of Wing 2 isten times that of Wing 1. 

In Fig. 7, the relation between flutter speed and A for Wing 2 has 
an apparent discontinuity near A = 15°. Increasing the value of GJ is 
found to have a pronounced effect on flutter speed at moderate sweep angles, 
but has little effect on flutter at high sweep angles. There is a drastic 
change in the flutter speed near 15 degrees sweep. For sweep angles beyond 
15°, the results obtained for this wing with high torsional stiffness resemble 
those obtained with the bending model and strip theory airloads. The reduced 
frequencies at the onset of flutter of Wing 2 are displayed in Table II for 
several sweep angles. From this table, it is seen that the flutter which 
occurs primarily as a body-freedom instability occurs at relatively low 
reduced frequenfcies when compared with the reduced frequencies which arise 
at the onset of bending-torsion flutter. Also, a comparison of the reduced 
frequencies in Tables I and II shows that the reduced frequencies at flutter 
in the two studies are comparable in magnitude. 
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Conclusions 


Before sumniarizing the results of this paper and listing conclusions, 
certain features of the idealized models studied should be reviewed. 

These models were chosen for analysis because of past experience with 
wind tunnel tests. A constant thickness, sheet-metal wing has a bending 
stiffness which is proportional to the wing chord measured perpendicular 
to the wing elastic axis; the torsional stiffness varies in a similar 
manner. This proportionality of the stiffness to the wing chord leads to 
bending and torsional stiffness distributions which are concave downward 
when plotted versus the spanwise coordinate. In actuality, the bending 
stiffness distribution which results from considerations of wing strength 
usually appears to have a concave upward distribution (cf. Ref. 3, p. 45). 

. The wings considered in this study had wing sectional centers of 
mass coincident, with the shear centers; thus, there was no elastic axis- 
c.g. offset.. Dynamic coupling was either nonexistent, as in the case of 
the uniform property wing, or minimal, as in the case of the variable 
property wing. This latter wing has a line of shear centers which is 
curved slightly forward when the wing is in its unswept position. 

The combinations and permutations of the various parameters which 
affect the aeroelastic stability of an aircraft are seemingly endless. 
However, several conclusions may be drawn from the present studies at 
zero Mach number. Prominent among these conclusions is that the inclu- 
sion of the rigid-body roll degree of freedom into the flutter model 
causes the critical mode of instability to change from an aperiodic in- 
stability (divergence) to an oscillatory instability (flutter). The 
degree to which the stability boundary is modified depends to a large 
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extent upon the sweep angle A and the ratio of the moments of Inertia 
In roll of the aircraft fuselage and the wing In its unswept position. 

If the 'wing instability appears as a coupling between wing bending- 
torsion deformation and rigid-body roll (the body->freedom mode), flutter 
speed is reduced as A increases. However, if the system parameters 
are such that flutter appears primarily as a fixed-root bending-torsion 
instability, the flutter speed may actually increase as the wing is sv/ept. 
If the wing can be either elastically or dynamically tailored, it may be 
possible to avoid the "body-freedom" mode type of instability altogether. 

Topics warranting further investigation include: the effect of Mach 
number on oblique wing, flutter; the significance of elastic axis - c.g. 
offset; and the effect of elastic tailoring of the wing. It is antici- 
pated that these and other studies will provide further insight into this 
unique aerodynamic design. 
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Table I - RcduKd Frequency k at Flutter (F1g. 2) 


k oic/2V^ 

A (Degrees)' 

0,0225 

15 

0,0330 

■ 30 

0,0494 

45 

0.0587 

■60 


Table II - Reduced Frequency k at Flutter 
(Wing 2, Fig. 7) 


k = WC/2V 

A (Degrees) 

0.29 

0. 

0.29 

7.5 

0.29 

15 

0.022 

20 

0.025 

25 

0.028 

30 

0.038 

45 • 

0.045 

60 
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Appendix 

The application of Galeirkin’s method to Eqs, (2,3) results in the 


determinant of Eq. 10. The elements of the matrix 


'll 


+ itiiD + 1 •!* g 

\%1 


d. . 

10 


are given below. 


(Al) 


d^2 = 0 


'13 


9 

8 


M- 


iwD 


(A2) 


(A3) 


dgi = 0 


^22 " “ + iwD + 1 - e 


\%) 


'23 


£ 

8 


y - ’ 


uD 


dg^ = w'^dSR/SO) - iw,|;(13/45) - (YSinA)(3/5) 


(A4) 

(A5) 

(A6) 

(A7) 


dgg " - /(13R/30) + i#(13/45) - (YSinA){3/5) 


(A8) 


ds3 = - u + iuij;(2/3) 


• (A9) 


The following definitions of terms are used in the above equations. 


D = (13/162-)(xL/VsinA) 
R 
3 


= (Iq/I-,^)cos‘^A 


= 5X/16 


tjj = yL/V 

3 

X = qcCj^^L sinAcosA/EI 
2 3 

Y = cos A/Ij 
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ELASTIC AXIS AND LINE OF 
AERODYNAMIC CENTERS 



Uniform Property King;' Definition of Geometrical Parameters. 




0 


Fig. 2 - Strip Theory Pre 
and Divergence S; 



QO - 163.5 MPH 
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Vc/V 



nn 


-Vd/V 


DO 


15 30 45 

A - DEGREES 


60 


Ratio of Aeroelastic Instability Velocity to Divergence Velocity at 
Zero Sweep Angle', Vj^^; Uniform Property Wing With Bending-Torsion 
Flexibility; Doublet-Lattice Aerodynamics. 




A - DEGREES 


Fig. 5 - Ratio of Aeroelastic Instability Velocity to Clamped Divergence 

Velocity at Eero Sweep Angle; Konuniform Wing Planform With Bending 
Torsion Flexibility. 
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Fig. 6 - The Effect on 



oqa/ja 



0 15 30 45 ^ 

A - DEGREES 


Fig. 7 - The Effect of Greatly Increasing the Torsional Stiffness Distribution 
of a Uniform Property Wing; Torsional Stiffness Distribution of Wing 

I 

2 is 10 Times That of Wing 1. 
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